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Abstract 

A model for describing structural pointlike defects in nanoscaled ferromagnetic materials is presented. Its de- 
tails are explicitly developed whenever interacting with a vortex-like state comprised in a thin nanodisk. Among 
others, our model yields results for the vortex equilibrium position under the influence of several defects along 
with an external magnetic field in good qualitative agreement with experiments. We also discuss how such defects 
may affect the vortex motion, like its gyrotropic oscillation and dynamical polarization reversal. 
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1 Introduction and Motivation 

Nanomagnetism has become one of the most promising branches in Nanoscience and Nanotechnology. Several 
nanostructured magnets have been fabricated and a great deal of efforts has been employed to investigate their 
physical properties, specially those believed to have potentiality for applications, for example, in magnetoelec- 
tronic mechanisms for data recording and storage, ultra-precise magnetic sensors, and so forth[TJ. 

Among several systems, nanomagnets whose shapes comprise vortex-like magnetization, as stable remanent 
states, have received considerable attention. Actually, such a pattern minimizes the total energy (exchange + 
magnetostatic) , for instance, in thin ferromagnetic nanodisks (with negligible anisotropy, like those made from 
Permalloy) [H [21 [H [5] . Several aspects of its structure and dynamics have been intensively investigated in the 
last years, like the interaction with litographically inserted defects [5J [71 [HI EE EH [HI Ell E4J> core polarization 
reversal [15l [TBI EO El, and so forth. For example, the presence of one or more (large enough) defects consid- 
erably modify vortex profile and dynamics; generally, its core is captured by one of the defects after some time, 
so that a pinned vortex appears to be the final state [6j El EH [HJ [12]. Core reversal, in turn, seems to be 
the one of the most important issues concerning vortex-like magnetization in nanomagnets, once the controlled 
reversion might open the doors for applications, like the utilization of the two polarization modes (up and down) 
as bits for recording and/or magnetic logic [TJ [T5l [T6J [17] . Nowadays, there are two distinct ways for achieving 
such a controlable reversal: applying suitable oscillating magnetic field pulses [15], or by means of spin polarized 
current [HI [17] . It has been additionally proposed, and observed in simulations, that polarization switching 
can be induced by vortex- hole interaction|13j or even by an interface dividing a nanodisk into two halves [14j . 
Actually, it has recently elucidated that core reversal takes place, at least in pure and homogeneous samples, 
whenever the vortex core achieves a critical velocity [18] , v cr ks 7 ^J2A/ no {A is the exchange stifness while 7 is 
the vortex gyromagnetic ratio; for typical Permalloy-made nanodisks A ~ 1CP 11 J/m and 7 ~ 10 4 Hzm/A, so 
that v cr ~ 10 2 m/s). Whether and how such a dynamical polarization switching is affected by both structural 
and litographically inserted defects remains untouched. 

Even though these and many other issues should be investigated for a better comprehension of such config- 
urations, it is noteworthy that even in the purest samples structural defects are present, generally randomly 
distributed throughout the material. In Ref. [19J . authors estimate in ~ 10 11 /cm 2 their density, in usual 
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Permalloy-made nanodisks. In this same work they also conclude that such defects may deeply modify vor- 
tex dynamics. In Ref. 20J, it has been carried out an experiment to analyze how those defects change vortex 
core path, as follows: An external field is suitably applied in order to displace the vortex core; if the sample 
were pure, a plot of the vortex equilibrium position against this applied field would give a perfect straight path. 
However, authors found a very complicated path with pronounced roughness along it. Such a roughness, together 
with successive jumps, were attributed to (small) structural defects. Here, we would like to present a proposal 
for modeling such defects as they were point-like. Although simple in its conception, our theoretical model yields 
a counterpart for the above-mentioned path which is in good qualitative agreement with experiments. We also 
discuss on the two possible types of defects, acting as pinning (attractive) or scatter (repulsive) site, and the 
important role they can play in both gyrotropic oscillations and dynamical reversion of polarization by means of 
their interaction with the vortex core. 



2 The model and its basic results 

The basic ingredient in modeling such defects is taking them as they were pointlike. Therefore, each defect is 
formally represented by a quantity like XiS 3 (x — x*i), where Xi is a real parameter, so that Xi > implies in a 
repulsive defect, while Xi < means that it is attractive. Varying Xi we may control the strength of the interaction 
with the defect located at xV If we focus on isotropic materials, then the total energy of a magnet is the sum of 
the exchange and magnetostatic contributions, as below: 

E nct = aJ (Vm) 2 dV + J rh{f) ■ , (1) 

where V is the volume of the magnet with normalized magnetization m(f) — M/M s = (m x , m y , m z ). Above 
the scalar demagnetizing potential (as usually, we are supposing the absence of free electric currents throughout 
the magnet) is given formally in terms of the effective magnetic charges, say, their volumetric, p m = V • M and 
superficial, a m = n ■ M, contributions, like below: 

* M s f V'-m(^) 1 . M s f n'-fhir 1 ),, 

47T jyi \r — r\ 47r J gi \r — r\ 

Of course, the explicitly expressions for both energies may be worked out if we fix the geometry of the mag- 
net and its magnetization configuration. In addition, magnetostatic contribution, coming from the long-range 
dipolar-dipolar interaction, is generally not easy to be evaluated. Here, we shall examplify our model by explicitly 
working out the influences of those pointlike defects whenever interacting with a vortex-type magnetization in a 
thin ferromagnetic nanodisk. In addition, our description is performed within the so-called rigid vortex regime, 
which is strictly valid if the vortex profile is not appreciably deformed, say, if its core is not largely displaced 
from the nanodisk center. Although simplifying considerably our analysis, it should be remarked that the main 
physical aspects of the relevant interaction will be brought about. 

Once the structural defects have been supposed to be pointlike, they will not appreciably alter the magnetostatic 
energy, since their volumes and areas, where effective magnetic charges should be computed, are strictly vanishing. 
This is further ensured provided that magnetization does not change abruptly around a given defect, which is the 
case for a vortex comprised in the magnet, as below. On the other hand, exchange term is modified, mainly at the 
vortex core, where such an energy is more concentrated. Actually, outside this core the defect-vortex interaction 
is approximately constant and much weaker than core-defect one. Thus, the important term in computing how 
the total vortex energy is modified whenever interacting with one or more defects is that concerning the core- 
defect interaction. Furthermore, once the vortex state will effectively interact only with those defects located 
on the nanodisk face, we essentially have a two-dimensional problem, so that A; is effectively measured in area 
units. Typical values for Xi parameter should take into account experimental facts, reproducing them at least 
qualitatively and, from the theoretical point of view Xi must be very small compared to the typical area unity, 
brought about by the exchange length, in a such a way that Xi « ig X . Such a relation is fulfilled in our analysis: 
while 1% X « 10~ 17 m 2 , usual values adopted here for Xi goes around 10~ 25 m 2 (see Figure 2 for further details). 
Then, while the defect-vortex interact, the parameter Xi show up indicating both its strength and if it is attractive 
or repulsive, as well. Eventually, such an interaction changes the exchange energy by AE ex = AL/2£j, with: 

6 = Qi{r- Si) I A^ 2 (x' - x i )(v 3 ,m(x') s j 2 d 2 x' . (3) 



2 



Here Qi is a function which accounts for how the defect located at xi interacts with the vortex core. Although 
its exact profile is difficult to be determined in practice, a gaussian function seems to be a reasonable choice, 
Gi(r~ £j) = e - a ( r -' x i) . [Other functions could be used to simulate such an interaction; disregarding its precise 
form, a free parameter, like a, should be introduced.] Its presence here essentially ensures that the vortex-defect 
interaction is smooth, starting effectively when the core border meets the defect, raising until a maximum, where 
the vortex center strikes the defect, then lowering. In addition, it plays also the role of controlling the interaction 
range to reasonable values, around the vortex core radius, Iq, in our case (see Figure 1). For Permalloy-made 
samples [3T], lo ~ 0.68/ OX (L/Z CX ) 1 / 3 w 8 — 15 nm, so that the range of this interaction S = 2lo/e ~ 10 nm, in 
agreement with experimental findings [19, 20]. Above, we have assumed disks with thickness L ~ 10 1 — 10 2 nm 
and l cx = \/2A/ hqM 2 ~ 5 — 6nm, once the saturation magnetization reads M s w8x 10 5 A/m, for this compound. 

In eq. ([3]), m = M/M s = (sin 9 cos </>; sin 9 sin cf>; cos 9) is the classical unit magnetization vector. In order 
to describe a vortex-type state with non-vanishing polarization, we must have <fi = 4>v = Q t&n~ 1 (y/x) + qn/2, 
where Q is the vorticity of the solution (topological winding number; Q is an integer) while q = ±1 is the vortex 
chirality (q = — 1 or g = +1 for clockwise or counter-clockwise circulation of the magnetization). On the other 
hand, (^-variable describing the vortex state is exactly known only asymptotically: cos9 v — > ±1 at the vortex 
center, r = 0, diminishing as one goes outside the core, where it vanishes. Among several suitable trial functions, 
and without loss of generality, we may take 9 V to be: 




if (4 » 

where p = ±1 is the vortex polarization (+ if the moment at the center points up or — if down). The parameter n 
may be choosen for adjusting the vortex core profile (see Ref.[55], for additional details). For simplicity, we shall 
take ?i = 1 in all of our calculations hereafter. Taking all these ingredients to eq. ^j, we readily obtain: 

f _ (2Zg-r')» + 4Za 

^~ Z4 (2Z g _ r2) (5) 

Therefore, the normalized energy density, u n ct = E net / AtthqMq V of the nanomagnet comprising a vortex inter- 
acting with a fixed defect at Xj, within the rigid vortex model, will read as below (the generalization for N such 
non-interacting defects may be readily get): 

u nct = |§ ln(l - s 2 ) [1 - &] + 2nF 1 (L/R) s 2 - h cxt s , (6) 

where R and L are the radius and thickness of the nanodisk. Parameter s measures the relative displacement of 
the vortex core from the geometrical center of the nanodisk, s = r/R (rigid vortex is ensured for small s; reasonale 
results still emerge at intermediary values, s ~ 0.2 — 0.3). Zeeman effect is accounted by the normalized field 
h ex t = H ext /M s . In turn, F\ is an auxiliary function, defined by: F\(z) = J °° J 2 (t) (e~ zt + zt — l)dt/zt 2 . Note 
also that the normalized energy is computed relative to the origin, s — 0, so that it does not give the precise value 
of the interaction between the vortex core and a defect located there. Here, such a value is normalized to zero, 
but it can be obtained explicitly by computing the exchange energy associated to the vortex interacting with a 
defect -Eex-int = E cx — E- mt = E eK (l — For our purposes, u nct is the relevant quantity. Now, expanding the 
exchange term in eq. ^ around s = and retaining only those contributions up to s 2 , we finally obtain that the 
equilibrium position of the vortex core, subject to an external field and interacting with a defect, reads: 

h 

(7) 



' cq 4irF 1 (L/R)-(l-&(l ex /R) 2 ' 
Clearly, expressions (JSC]) recover their usual forms [S] as long as we remove the defect, = 0. 



3 Further results and discussion 

Now, let us consider a number of defects, for instance along a given line, and apply a suitable field for moving the 
vortex core towards each defect. This will enables us to show how our model describes vortex-defect interaction 
more clearly, making possible a comparison with available experimental findings. For that, let us place several 
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pointlike fixed defects, with different interaction parameters A, along x-axis. An external field is then applied 
to move vortex core along this line, so that it interacts with these defects while moving. The resulting plot of 
its equilibrium position against the applied field is shown in Figure 2. Note that, by virtue of the interaction 
with the defects, the vortex core trajectory is no longer a straight line, which is the path for a pure sample, with 
£i = in eq. ||7J). Besides of being quite complicated it should be noted that this actual trajectory presents 
a remarkable roughness, mainly around each defect position. We also note two distinct sorts of jumps in this 
plot: i) abrupt changes in s cq , at practically the same value of h cxt , what says us that the core is interact- 
ing with a repulsive defect, A > (indicated by black arrows in Fig. 2); ii) conversely, in other places s oq is 
kept practically unaltered while h cxt changes considerably, evidencing the interaction with an attractive defect, 
A < (as blue arrows indicate in Fig. 2). Comparing these with the experimental results provided in Ref. [20] . 
mainly with Figure 5 and related text from this article, a good qualitative agreement between them can be realized. 

Once the vortex core position is changed whenever interacting with scatter or pinning defects throughout the 
sample, we may wonder whether and how these structural defects could modify the gyrotropic motion and its 
associated frequency. Qualitatively, it is expected that its traced path will be much more complicated than that 
smooth ones in pure samples. The roughness along its actual path will resembles us a Brownian-type motion, 
enlarging the distance followed by the vortex to complete a revolution. If the concentration of pinning and scatter 
sites is similar each other, a net compensation is expected so that the gyrotropic frequency is expected to remain 
practically unaffected, provided that the scattering does not enlarge the trajectory of the core. However, if we 
have much more, say, pinning sites, a net retardation in the vortex motion should be verified causing a decreasing 
in the frequency. Such an effect could be estimated qualitatively, for instance, if we consider the Landau-Lifshitz 
equation: 

G x v = kr + m v a , 

with the reasonable assumption that the vortex mass is very small [12] , m v w 0. In this equation, G = 
— (Mg/7) J v V(cos0„) x V<fi v dV is the gyrotropic vector of the vortex, with 6 V and <j> v being the vortex state 
functions, previously given. In the simplest case of a vortex in a circular nanodisk it may be evaluated to give 
G = (2tt M s L/j)pqz, evidencing its topological nature. Parameters r, v and a are the position, velocity and 
acceleration of the vortex core, while k plays the role of a spring-like constant, derived from the net energy of 
the magnet as usual, F = —kr = —VE nct . In this case, supposing small vortex core displacement, s << 1, k is 
explicitly obtained to be: 

I 2 



k = ttLM. 



A-kF\{L/ R) — (1 - £i)^p 



Once k is changed by the ^-factor, we clearly see that the structural defects are expected to modify the gyrotropic 
motion of the vortex core. However, we cannot proceed further with our analysis because we can no longer 
assume that the vortex core perform harmonic oscillations around its equilibrium position (as is usually done), 
once deviations may be appreciable whenever defects are present, making our analytical approach efficientless; by 
virtue of that, numerical simulation appears to be a more suitable way for studying how these defects affect the 
gyrotropic frequency [23] ■ Moreover, they can also change the situation for core reversal take place. Actually, 
even assuming that defects does not change the critical velocity (in pure samples, it only depends on the exchange 
stiffness, v cr ~ r yy/2A/no [18]) they may considerably change the scenario of polarization switching: For instance, 
if vortex core is moving in a region with larger concentration of attractive defects, the successive pinnings will 
deaccelerate its core, dissipating kinetic energy as spin waves and/or heat, preventing its immediate achieving of 
the critical velocity. Similarly to the case of gyrotropic motion, a precise answer demands further investigation, 
namely, by means of numerical simulations. 



4 Conclusions and Prospects 

A model for structural defects in nanoscaled ferromagnetic systems is presented. Once such defects are described 
as pointlike objetcs their main modification appears in the exchange contribution to the total energy. As an ex- 
ample, we have worked out the case of a vortex-type pattern comprised in a thin nanodisk. An expression for its 
equilibrium position whenever subject to both, an external field and defects, is analytically obtained, immediately 
recovering its usual counterpart if the defects influence is removed. A good qualitative agreement between experi- 
mental and our theoretical results, concerning how structural defects modify the vortex motion, is clearly realized. 
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Among other prospects for future investigation, we may quote the refinement of some ingredients entering this 
model (defect-magnetization interaction form, strength of the interaction, and so forth), for instance, to explicitly 
work out how such defects affect other magnetization configuration and/or sample geometries. We also intend to 
reconsider this model, but for studying magnetocristaline anisotropy. In this case, we expect that instead of A.; we 
should consider a vector function, A(af), which should account for the preference of the magnetization in pointing 
in a given direction at each point or region inside the nanomagnet. Such a study could perhaps shed some extra 
light on the pronounced asymmetry in the gyrotropic frequency between up and down polarizations, as observed 
in the experiments of Ref . [H] , and very recently addressed in the work of Ref. [55] . 
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Figure Captions 

Figure 1: Illustrates how the defect-vortex interaction takes place: tiny at the core borders, where exchange 
energy is small, raising to its maximum at the vortex center, where exchange density gets its highest value. 
At r = the interaction is maximum whose strength reads AXi/Iq. Its effective range S is around 10 nm, in 
accordance with recent experiments; its (Gaussian) shape may be controlled by varying a-parameter (we 
have taken oT 1 = 2.5 x 10 4 in order to set the interaction range ~ 10 nm, according to experiments). 

Figure 2 (Color online): How the vortex equilibrium position versus external field is modified by virtue of the 
interaction with structural defects (if the sample were pure, a straight line would be the case). A comparison 
with experimental findings, as reported in Ref. [20] , shows a good qualitative agreement. Namely, note the 
interaction with pinning and repulsive defects, indicated by blue and black arrows, respectively. [Typical 
values used for |Ai| read in the range ~ 10~ 25 — 10~ 26 m 2 . Note that such values are much smaller than 
^ 2 X ~ 10~ 17 m 2 , the elementary exchange area, showing the smallness of our described defects, as initially 
supposed] . 
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